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This paper addresses the question how to implement a desired two-qubit gate U using a given 
tunable two-qubit entangling interaction Ti.i„t(t). We present a general method which is based 
on the K\AKi decomposition of unitary matrices G S'C/(4) to calculate the smallest number of 
two-qubit gates Ui n t(t) [based on Tiint(t)] and single-qubit rotations, and the explicit sequence of 
these operations that are required to implement U. We illustrate our protocol by calculating the 
implementation of (1) the transformation from standard basis to Bell basis, (2) the CNOT gate, and 
(3) the quantum Fourier transform for two kinds of interaction - Heisenberg exchange interaction 
and quantum inductive coupling - and discuss the relevance of our results for solid-state qubits. 
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I. INTRODUCTION 

How to implement two-qubit quantum circuits in an optimal way, by which we mean here that they require a 
minimum amount of operations, has become a topic of active interest over the last years. Typically, the central 
question asked is how to perform an arbitrary operation £ SU (4) with the least number of two-qubit gates such as 
the CNOT (controlled NOT) gate. This has led to the discovery that any two-qubit operation can be implemented 
using three CNOT gates and local gates in between^. Alternatively, two-qubit operations can be constructed by three 
applications of the (SWAP)" gate^, with < a < 1, or by two applications of the B-gate^. From a practical point of 
view, it is essential to perform quantum operations as fast as possible, before decoherence (due to interactions with 
the environment) sets in. This is of particular importance for solid-state qubits, which in general interact strongly 
with their environment. For a specific physical system, however, efficient implementation of quantum operations does 
not necessarily coincide with minimizing the number of CNOT gates, if CNOT cannot be directly generated by the 
given two-qubit Hamiltonian. Rather, efficient construction of quantum operations requires minimizing the number 
of two-qubit gates which are based on the (entangling) two-qubit interaction that is naturally available in the system, 
as well as minimizing the number of single qubit operations (which in practice generally also take a finite amount of 
time to implement^). 

In this paper, we present a general, explicit and self-contained analytical protocol for calculating the implementation 
of a desired two-qubit quantum gate U £ SU(4) as a sequence of single-qubit rotations and two-qubit gates Ui nt {t), 
where Umtit) is based on the two-qubit interaction H.i n t(t) that is naturally available in a given system, Ui n t(t) = 
exp[— ^ J dt Hi n t(t )]■ We assume that U int {t) has some entangling capacity, which is a necessary requirement if U 
itself is an entangling gate. We also assume that Umt(t) is a tunable gate, i.e. that one can switch the interaction 
'Hintit) in the system on and off during well-defined amounts of time by changing external parameters. This is often 
the case in real qubit systems where the qubits consist of massive particles such as electrons or ions. Our scheme is 
based on separating the two-qubit (entangling) part of U from the single-qubit basis transformations (the so-called 
K\AK 2 or single-value decomposition^, with A G SU(4) and Ki,K 2 € SU(2) x SU(2)), translating each of the 
matrices K\, A, and K 2 into a sequence of operations which involves the smallest possible number of U% n t{t) for a 
suitably chosen time t, and minimizing the number of single-qubit rotations in the total sequence of operations using 
permutation relations. The resulting shortest-sequence implementation of U - which we refer in this paper to as the 
optimal implementation ofU- contains the least number of operations from the given " library" of single- and two-qubit 
gates 7 . We illustrate our scheme by calculating the optimal implementation of three elementary quantum operations 
U [the transformation from standard basis to Bell basis, the CNOT gate, and the quantum Fourier transform] for 
two kinds of two-qubit interaction Hint that are relevant for solid-state qubits: Heisenberg exchange interaction 
(corresponding to the (SWAP) Q -gate) and quantum inductive coupling (corresponding to the B-gate). Either of 
these gates, in combination with single-qubit rotations, forms a so-called universal set^^i^ into which any operator 
U € SU (4) can be decomposed. 

The paper is organized as follows. In Sec. [IT] we outline the decomposition protocol, starting with the decomposition 
of an arbitrary gate U € SU(4) into matrices K\, A and K 2 . Sec. Ill Bl addresses the special case in which U is a 
maximally entangling gate, for which the protocol can be simplified. In Sec. Ill CI we show how each of the matrices 
K\ , A and K 2 can be decomposed into a sequence of operations consisting of single-qubit rotations and a minimum 
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number of two-qubit interactions Ui n t(t) of a given kind. In Sec. Ill Dl we show how to iteratively minimize the 
total number of single-qubit rotations in the decomposition of U, and thereby obtain the implementation of U that 
involves the smallest number of both Ui n t(t) and single-qubit rotations. In section [Hi) we illustrate the decomposition 
protocol for two kinds of interaction that occur in solid-state qubit systems: Heisenberg exchange interaction Hex 
and quantum inductive coupling Timd- We show how to implement the transformation from standard basis to Bell 
basis (Sec. IIII Aj) . the CNOT-gate (Sec. IIIIB[) and the quantum Fourier transform (Sec. IIII Cj) using either of these 
interactions. Finally, in Sec. IIVI we discuss how our results can be implemented for electron spin qubits in quantum 
dots and for superconducting flux qubits, followed by conclusions in Sec. fVl 



II. EXPLICIT DECOMPOSITION PROTOCOL 



A. General gates U £ SU(4) 

The decomposition protocol we present in this section is based on the K1AK2 decomposition^, which states that 
every two-qubit operation U £ SU(4) can be written as U = K1AK2, with Ki, K 2 E SU(2) (8 SU(2) and 



A = exp [i(aa x <g> a x + (3a y ® a y + <g> a z )\ 
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with = cos(a ± (3) and Q^g = sin(a ± 0). Here a, (3, 7 € [0, 27r) and <j x , a y , a z denote the Pauli matrices. Our 
goals in this section are 1. to calculate the matrices K\, K% and A for a given U £ SU(4) and 2. to translate each of 
these matrices into the smallest number of single-qubit rotations and two-qubit interactions Ui n t(t) of a given kind. 
In this subsection and the next we develop a systematic protocol for the first goal 9 . The implementation of Ki, A 
and K2 using Umt{t) and single-qubit rotations is the subject of section Hi CI 
We begin by introducing the matrix Q G SU (4) defined as 



10 i 

l_ I i 1 

V2 i -1 

1 -i 



(2) 



Q transforms the standard basis { |ei) , |e 2 ), |es), |e4)} of SU(4) into the magic basis {|mi), \m 2 ), \ms), \1ri4,)}. Here 
the standard basis is defined in spin-1/2 notation with the spin along the z-direction as \e\) = |tt)j l e 2) = |TI), 
|e3) = ||t): l e 4) = III) an d the magic basis is defined as: 



\mi) 
\m 2 ) 
|m 3 ) 
|m 4 ) 



1 

1 
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(ITT) + III)) 
KIT!) + lit)) 



^(ITD-IIT)) 
^jf (ITT) -III)) 



(3a) 
(3b) 
(3c) 
(3d) 



Next, we define 



Ox = Q^KxQ 
2 = Q^K 2 Q 
F = Q^AQ, 



(4) 
(5) 
(6) 
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where Ox, 2 G SO (4)12. and F = diag(Ai, A 2 , A 3 , A 4 ) with \ u . . . A 4 G C. Using Q, © and ©, we find that the 
relation between Ai, . . . , A4 and a, 8, 7 is: 

\ 1 = e *(«-/J+7) ( 7a ) 

A 2 = e^ *" - ^ (7b) 
A 3 = e ^- a -P-^) (7c) 
A 4 = e 4( - Q+/3+7) . (7d) 

Our aimii is to find the matrices 0\, O2 and F, and from those via (HJ)-© the matrices K\, K2 and A. To this end, 
we proceed in 3 steps: 

Step 1: We first calculate Umb, where Umb is the representation of U in the magic basis ©: 

U MB = Q f UQ = OxF0 2 . (8) 

We then define vectors \vi) as 

|w 4 > = O^lei) fori=1...4. (9) 

The vectors (O are also the eigenvectors of U mb Umb with corresponding eigenvalues fXi = Af : 

Ul IB U MB \vi) = T 2 F 2 2 \ Vi ) = OlF 2 \e % ) = \ 2 O t 2 |e 4 ) = W . 

Thus by calculating the eigenvalues and eigenvectors of U mb Umb, which is a known matrix, we directly obtain the 
matrix O2 from Eq. © as 

2 = (h)..>4» T 

and hence, using Eq. ([5]), 

K 2 = Q(\v 1 )...\v 4 )) T Ql (10) 

In the above we are free to assign the eigenvectors of U mb Umb m some chosen order to the vectors \vi), with the 
restriction that since O2 G SO (4) it should hold that detC>2 = 1. 

Step 2: Using the eigenvalues Af of U mb Umb found in Step 1, we calculate the arguments a, (3 and 7 of 
the matrix A [Eq. §\$\ 12 using £[]) : 

—i 1 
a = —\og(\ 1 \ 2 \- 1 \- 1 ) + -k 1 TT (11a) 

3 = Z!l g(A 1 - 1 A2A3 1 A 4 ) + ifc27r (lib) 

7 = log (AiA 2 _1 A3 X A 4 ) + 7^371") (11c) 

where k\, k 2 , k% G Z. 

Step 3: Finally, the matrix 0\ is calculated from 



Oi - Oi^0 2 OjF* = UmbOIf", 



and hence 



#1 = QUmbOIF*Q? (12a) 
= U(QQ T )K 2 A*(QQ T )*. (12b) 

There are two degrees of freedom in this decomposition procedure. First, there are different ways of identifying 
the eigenvectors of U mb Umb with the vectors \vi). Secondly, there is some freedom in choosing the arguments a, 
(3 and 7 of the operation A. In general, different choices lead to different matrices Ki, A and K 2 , and the KiAK 2 - 
decomposition of U is not unique. 
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At this point it is useful to compare the procedure for obtaining a .FLiAi^-decomposition described above to the 
one proposed by Kraus and Cirac (KC) 5 . There are two main differences between the two: 1) Using KC's method 
the matrix K2 is given in the form SU(2) x SU(2) and involves 3 unknowns (Euler angles), and 2) KC's method 
leads to fixed values of the Aj's (instead of freedom of choice of the sign of each Aj). Which of the two methods is 
more convenient to use depends on the goal of the calculation to be made. If this goal is to (only) find a K1AK2 
decomposition of a given matrix U, then having if 2 in SU(2) x SU(2) form and the A;'s fixed may be advantageous, 
if it leads to less equations to be solved. If, on the other hand, the goal of the calculation is to first find a K\AK~i 
decomposition of U and then decompose each of the matrices K\, A and Ki in terms of single-qubit rotations and 
a given tunable two-qubit interaction Ui n t(t) (as it is in this paper), the freedom of assigning the vectors \vi) and 
choosing the sign of the Aj's that our method entails is advantageous, since it allows us to choose these parameters in 
such a way that the ensuing decomposition into single-qubit rotations and Ui nt (t) is straightforward. 

B. Maximally entangling gates U £ SU(4) 

The procedure described in the previous subsection can be simplified if U is a quantum gate which maximally 
entangles a suitably chosen separable basis, which we refer to as a "maximally entangling operator" or a "quantum 
gate with maximally entangling capacity" (not to be confused with quantum gates which "only" maximally entangle 
a single separable quantum state instead of a whole basis). For this special class of two-qubit operations the following 
Proposition holds: 

Proposition 1: Any transformation M € SU(4) from the standard basis to a maximally-entangled basis can be 
expressed in the following way: 

M = K t e X«r-®«r. + «r,0«r,) ^(1) (^(2) {ly)R (l) ^ R (2) (^(1) ^(2) ^ {u) 

where e, fi, is, £, £, 77, and u> satisfy one of the following three conditions: 

1) e = 0, C,C€{0,±7r}, ftI /,t),we [0,2tt) (14a) 

2) < e < 1, fi, u, £, C € {0, ±tt}, 77, lu £ [0, 2tt) (14b) 

3) e=l, M,f€{0,±7r}, f, C, V, " e [0, 2tt). (14c) 

Proof: It has recently been shown 13 that any maximally entangling operator is of the form A(e) = 
exp [(l/4)in(cT x ® cr x + ea z cr z )] (with e £ [0,1] and modulo permutations in o~ x , o~ y , o~ z ). We are thus left 
with the question which rotations may be performed on the standard basis before applying A(e), such that A(e) 
acts as a maximally entangling operation. Since rotations around the z-axis only multiply each quantum state by a 
phasefactor, they are always allowed. It is then necessary and sufficient to investigate for which values of e, /1, v, £ and 
£ the transformation T = A(t)R x 1 \pL)R x 2 \v)R^y\^)R l 'y\( > ) corresponds to a maximally entangling operation. These 
values can be found in a straightforward way by determining for which values of e, /i, v, £ and C, the concurrence^ 
equals 1 for each of the four column vectors \wi). The calculation is given in the Appendix and leads to the solution 
Eq. (HID. ■ 

Note that in order to transform from the standard basis to a specific maximally entangled basis the conditions on 
the angles in Eq. (|14[) are in general more restrictive. For example, to transform to the Bell basis we require /i, v, £, £ 
£ {0,±7r}, 77, uj 6 [0,27r) and cti,f3i,ji € [0, ±7r/2, ±tt] for all values of e. 

For gates U £ SU(4) which maximally entangle the standard basis, step 2 in the decomposition procedure 
from the previous subsection can now be simplified: using Proposition 1 we choose A1-A4 such that j3 = and 
a = 7r/4. For each of the possible assignments we then calculate 7 (and thus the value of e in Eq. (jT3j) ) from 
Eq. (TTTcTl . 

Having found the decomposition of U into the matrices K\, K% and A, we now proceed to develop a systematic 
scheme to calculate the implementation of each of these matrices using a sequence of single-qubit rotations (for K\ 
and K2) plus two-qubit interactions Ui n t(t) of a given kind (for A). This is the topic of the next subsection. 

C. Decomposing Ki, K% and A into Ui„t(t) and single-qubit rotations 

We first consider gates U £ SU (4) with maximally entangling capacity, as defined in the previous subsection. 
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Decomposition of K% into single-qubit rotations. - For maximally entangling two-qubit gates, a decomposition 
of the matrix K 2 into single-qubit rotations can be found in a straightforward way by comparing the matrix K 2 
obtained from (fit)]) to the sequence of rotations on the right-hand side of Eq. (fl3|) , using the appropriate conditions 
on the angles (jT5J) that corresponds to the value of e = A'y/ir calculated from Eq. (|llc[) . The resulting decomposition 
of K 2 is not unique, as can be seen from (| 14[) where different sequences of single-qubit rotations each correspond to 
a correct decomposition of Ki- 

Decomposition of K\ into single-qubit rotations. - K\ can be translated into single-qubit rotations using the 
Euler decomposition theorem-^, which states that every K S SU{2) ® SU(2) can be written as 

K = [] R ( J ] 4 3) (Pi) R i j) (7i)- (15) 

Here Rz(aj) represents a rotation of qubit j over ay around the z-axis and Oj, f3j, jj € [0, 2tt). By substituting i^i 
calculated from Eq. (TT!?|) into the left-hand side of Eq. (fT5|) and solving the resulting set of equations, we directly 
obtain the values of the angles Oi, (3i and 7$. It is clear that, as for K2, the decompositon of K\ is not unique and 
that different sequences of rotations all yield a correct decomposition of K\ . 

Decomposition of A into Ui„t(t) and single-qubit rotations. - Our next goal is to decompose the matrix A 
[Eq. ([T]) with the values of a, (3 and 7 obtained from Eq. (TTTj) ] into a sequence of two-qubit interactions Ui nt (t) and 
single-qubit rotations such that the smallest number of Ui n t(t) is used. To this end, we first examine the entangling 
capacity of the given two-qubit gate Ui n t(t) as a function of t, i.e. we investigate how many applications of Ui n t{t) 
are needed to transform the standard basis into a maximally entangled basis for a given t. For tunable two-qubit 
gates Ui n t{t) there exists in general a value of t, say t* , for which this minimum number of applications is small. 
The value of t* is obtained by calculating a measure of bipartite entanglement, such as the concurrence^, for Ui nt {t) 
applied to the standard basis vectors. In many cases t* can be determined directly from inspection of Ui n t(t) (for 
Ui n t{t) written in the standard basis). For example, if Ui nt {t) is the (SWAP) Q -gate (see Sec. Mil Eq. ([T7|) ). t* is 
determined by a. In this case two applications of (SWAP) 1 / 2 (plus a single-qubit rotation in between) are sufficient 
to implement ? 7 16 ' 17 , since (SWAP) 1 / 2 maximally entangles two out of four standard basis states (and affects the 
other two basis states only trivially, by multiplication with a phase factor). Therefore t* = ir/(2h~J) (for constant J, 
see Eq. (|16l) below). 

Once a sequence of operations, say A, that transforms the standard basis into a maximally entangled basis and 
contains the smallest possible number of Ui n t (plus single-qubit rotations) has been found, A and A differ by local 
operations only. The latter can be calculated in the same way as for K\, by using the Euler decomposition (fT5|) of 
AA~ l . 

For two-qubit entangling gates U € SU(A) that do not have maximally entangling capacity, the decomposition 
procedure of K±, A and K2 into Ui n t(t) and single-qubit rotations is the same as for maximally entangling gates 
described above - except that the matrix K2 now also has to be decomposed by using the Euler decomposition (|15[) . 
since Proposition 1 does not apply. 



D. Optimization 

By joining together the individual decompositions of K\, A and K2 obtained via the method described in the 
previous subsection, we have obtained a decomposition of U = K\AK-x which consists of the smallest number of 
two-qubit interactions (plus a number of single-qubit rotations). This optimal decomposition is obtained by first 
optimizing the number of two-qubit gates and subsequently optimizing the number of single-qubit gates. The total 
number of rotations in this decomposition can often be reduced by using one or both of the following: 

(i) Commutation relations between Ui n t(t) (or A, as defined in the previous subsection) and rotations, 

(ii) Euler's theorem: writing rotations Rml(4>) asi£ Rn}{a\) Rp} '(</>) Rn} {a.2) with m 7^ Pi, both m and pi orthogonal 
to rrii, and 0.1,0.2 € {±7r/2} (and using (i) again). 

Once both of these do not lead to a further reduction of the number of single-qubit rotations, the minimum number 
of rotations has been found. 

In practice, depending on experimental conditions, it may only be possible (or be easier) to implement rotations 
around certain axes and therefore be necessary to e.g. translate rotations around the z-axis into a sequence of 
rotations around axes in the (x, y)-plane. By implementing specific requirements such as these and using commutation 
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relations to minimize the number of rotations we obtain the implementation of U which requires the smallest number 
of operations from a given library of single- and two-qubit operations. 



III. ILLUSTRATION OF THE DECOMPOSITION PROTOCOL 



In this section we illustrate the protocol developed in the previous section for two types of interaction Ui nt {t): 
Heisenberg exchange interaction and quantum inductive coupling. Heisenberg exchange interaction is described by 
the Hamiltonian 

n EX (t) = (i/A)h 2 j(t)d^ -<j^\ 



with J(t) the time-dependent (tunable) exchange energy. 
(SWAP)" = U E x(t) = exp[-(i/fi)/ * H B x(r)dr] with 



This interaction corresponds to the (SWAP)"-gate: 



a(t) = / J(r)di 

n Jo 



or, equivalcntly, 



(16) 



(SWAP)" = e 



e% i7r \ 

cos(|tt) isin(fTr) 
zsin(fTr) cos(§tt) 



V 



o 



o 



(17) 



/ 



The second type of interaction we use to illustrate our protocol is a tunable coupling between two magnetic fluxes 
(see Section IIVI for a description of a physical realization using so-called flux qubits) . The Hamiltonian corresponding 
to this type of interaction is given byi^ 



rtind(t) — — [cr x 



o y ) — (a x 



><7y), 



(18) 



7±(i) = \ Jq 9±{T)dT, where g±(t) are tunable system parameters. The Hamiltonian (|18p corresponds to the so-called 
B-gate^: £(7+,7_) = U int {t) = exp[--| f Q H in t(r)dr] or, equivalently, 



'7+ 



■S(7+,7-) = 



/ cos ( 2 > 




Visin(^) 





cos(^-) 
i sin(^-) 






i sin(^-) 
cos(i-) 



zsin(^j-) 



COS(^L) 



(19) 



with 7± G [0, 27r). Note that Eq. (fTi?)) for 7+,7- = ±7r/2 maximally entangles the entire standard basis, whereas the 
exchange interaction \/SWAP only produces maximal entanglement for two out of the four states in the standard basis. 
For both H,Ex{t) and 7ii n d(i) we now calculate the optimal (shortest-sequence) implementation of three elementary 
quantum operations: the transformation from standard basis to Bell basis, the CNOT gate, and the quantum Fourier 
transform. The former two are maximally entangling operations, but the latter, as shown below, is not. 



A. Transformation from standard to Bell basis 



A transformation from the standard basis to the Bell basis is by definition a maximally entangling operation. It has 
been shown that when using the Heisenberg exchange interaction [The (SWAP) Q -gate, Eq. (fTT]) ] the shortest sequence 
of operations that transforms the standard basis into a Bell basis is given by^^ 

M$ = VSWAP (tt) VSWAP 



/ 









A 


e # 




















e -i<t> 


V 





-ie lrt > 


e i4> 


o / 



(20) 
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where j=l,2 labels the qubit and B^} (jr) represents a rotation of qubit j around an arbitrary axis = (cos </>, sin <j>, 0) 
in the (x, y)-plane. The analogous transformation for the quantum inductive coupling consists of a single application 20 
of the i?-gate [Eq. (fll?)) ] with 7+, 7_ = ±S. In order to find the transformation from the standard basis to the 
"standard Bell basis" {(l/\/2)(|TT) ± I U», (l/\/2)(lti> ± I IT))}, we need to decompose the matrix 



If 



st->Boll 



1 
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/ 1 o o 1 
Oil o 
1-10 

V i o o -l 



(21) 



(or a permutation of (|21[) in which the columns are interchanged). Using the protocol outlined in the previous 
section, we first calculate the matrices ifi, A, and if 2 and then decompose these matrices into the shortest sequence 
of Ui n t — (SWAP) Q or Uint = B plus single-qubit rotations. 

Starting with Step 1 in Section III Al we find that 



st— KBell 
MB 




and the eigenvalues and corresponding eigenvectors of (U M ~^ ^mb ) are gi ven by 

Mi = Ms = i, 

(i 3 = M4 = -i, 1 03) 
We now choose |uj) = |oi)Vi=l,...,4 and calculate the matrix K 2 using ((TI 




(22) 



if: 



st-»Bcll 












l + i 








-1 -i 
















1 - i 












(23) 



Next, we use the eigenvalues Aj = i = 1, ... ,4 [Eq. (f22l ] to calculate the matrix A from Eqns. (fTTj) and |T]). 

Choosing Ai = v?j A2 = — yi, A3 = >/— ?, A4 = —\J~-l we find the solution a = 7r/4, /3 = 7 = and hence 



A 



st->BcH 




(24) 



Finally, we obtain K\ from Eq. (fT2|) : 



if 



st->BeU 




(25) 



We now decompose each of the matrices ifif^ 3611 , yi st ^ Bcn an( \ Kf~* Bo11 [Eqns 

qubit rotations and (SWAP) Q -gates, using the procedure described in Section III CI Starting with ^st-+Bcii anc j 
comparing the matrix |23|) to the right-hand side of Eq. (|13|) for condition ()14a|) . we find 



into a sequence of single- 
if 



(26) 
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By comparing the matrices in Eqns. (|24|) and (f20|) . we directly obtain a decomposition of v4. st— >Bell ) since for — x 
and j = 1 the two matrices only differ by a spin flip of the first qubit (which exchanges the first and third row as well 
as the second and fourth row). Hence we find, disregarding a global phasefactor, 

^st^Bdi = R (i)^) VSWAPi?«(7r) VSWAP. (27) 

This decomposition of ^4 st ^ Bo11 [ s a \ so obtained by following the decomposition procedure described in Sec. Ill C[ 
by first calculating the entangling capacity of (SWAP)" Va. Finally, the decomposition of the matrix Kf 1 ^ E 
SU(2) ® ST/ (2) is obtained by using Eq. fTS]), and we find: 

#f-Beii = ^(i) ( _ZL )i? (2) (7r) _ (28) 



The total decomposition of [/ st -> Bc11 i n terms of V S WAP-operations then becomes: 



jjst-tBcll _ ^st^Bell^jst-^Bcll^sWBcll 

= R { P (™) W i£> W ^/SWAP i£) (tt) VSWAP i?« (tt) i?f (tt) J#> (|) 

= ^H-^VS^i?^^)^^^^ 1 )^)^ 1 ^!) (29) 

In the last step of Eq. (J2SJ) we have used the relation VSWAP (tt) VSWAP uf 1 (tt) = 

R ( z ] (tt) (tt) VSWAP R { } ] (tt) VSWAP. Eq. (29) is tnc shortest sequence of operations that can be used to 
implement [/ st ^ B ° n using Heisenberg exchange interaction and single-qubit rotations. Repeating the same procedure 
for other transformations from the standard basis to the Bell basis (obtained by permutations of the columns of 
Eq. we find the following possible optimal decompositions of [/ st - >Boll : 

rrrt^Bell _ J {~<J f )R^ ] i#> (± § ) 

^(swap). - i fl (i) ((jf ) 4D (7r)M W jR W (±f) ^ 



i?i 1) (-af)Mi i) i?^( 7 r)i?i 1) (±2, 



with (i, a) G {(1,1), (2, -1)^1. 

The analogue of Eq. ([50)) for quantum induced coupling (the B-gate, Eq. §T§§) is given by: 

' R£\-a%)B(a%,aZ)R£\±iL) 
M 1) ( ( r|)i?i 1) WB( ( 7§, ( Tf)i?i 1) (±f) 

i^ 1) (-<7f)iii 1) (7r)B((r|,£rf)i^i(7r)M 1) (±f) (31) 
fli 13 (<xf W5(<t| , df )i?i X) (vr)i?i 1) (±§) 
i?i 1) ( ( 7§)i?(af,af)i?i 1) (7r)i?i 1) (±f) 

with a = ±1. Note that Eqns. (|50)) and (|3 1 [) involve six rotations more than the decomposition of M n ,, Eq. (|20[) . 
From a practical point of view, it is thus more efficient to implement a transformation from the standard basis to a 
Bell basis with complex coefficients than to the "standard Bell basis" when either Ui n t =(SWAP) Q or Umt = B is 
used. 

To conclude this section, we also give the decomposition of two alternative matrices ^4 st ^ Bc11 . The first one is 
obtained by using the K\AKi- decomposition procedure of Kraus and Cirac^ and reads: 




^ait7 - | n _,■ n n I ■ ( 32 ) 
^aJti B ° n does not maximally entangle the standard basis. In order to find its decomposition, wc first calculate an 



unentangled basis for which ^jfc* does act as a maximal entangler. To this end, we note from the matrix 
that the latter basis must be obtained from the standard basis by creating a superposition of both the first and the 
second qubit, since Eq. couples | ff) to | J.J.) and | fj.) to | \\), Mathematically, a possible choice of rotations is 



1-1-11 
1 1-1-1 
-11-11 
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which translates into (using Euler's decomposition theorem) R^\j) R* (it) Ry*\— -|). Thus A| t 1 jTj >Be11 

R ( y ] '(f) R ( z ' '(tt) i?^' '(-§) = (basis rotation) x VSWAP (tt) VSWAP ee (basis rotation) x m| 1} . The remain- 
ing basis rotation is easily found by using Euler's theorem and the full decomposition then reads 

^u Bc11 = b$X\ ) 4 2) (f ) M P *P 00 4 2) (f )• 

As a second example of an alternative matrix A^*t2 B > we consider the matrix obtained by using a different choice 
of Aj's, namely Ai = X2 = — V*, A3 = yi and A4 = — >/— ' 1. For this choice the matrix A [Eq. (JTJ] becomes 

/ 1 \ 

4 st->Boii _ I -i 1 r . 

Aat2 - o o _i ' w 
\o 1/ 

^ait2 B ° n ac ^ s as a maximal entangler on the same unentangled basis as .Ai^j* Bo11 (see above), so that A s ^ Bcl1 
R { y ] {^) R { z ] {tt) R { y ] {-^) = (basis rotation) x Af| 1} . The remaining basis rotation is again found by decompos- 
ing Af^ BeU i?^ 1) (f)i?i 1) (7r)i?J) 2) (-f) (M^)- 1 using Euler's theorem and we find A^ Be11 = ^(-f) Ry\-%) 

Mi 1] b£\*) 4 1} (-d 4 2) (D- 



B. The CNOT gate 

The CNOT-gate is given by 

/ 1 
CNOT<"' _ ° I I 1 
\ Q 1 . 

In order to find the optimal (shortest-sequence) decomposition of the CNOT gate in terms of the (SWAP)" and 
B-gate, we again use the protocol from Sec. Inland proceed in the same way as for [/ st -* B ° n i n the previous subsection. 
We first calculate K\, A and K2, then decompose each of these matrices into (SWAP) Q or B-operations plus 
single-qubit rotations, and subsequently use permutation relations to optimize the number of rotations. 

Starting with Step 1 in Section Hi A[ we find for the representation of CNOT 1 - 1 ' 2 '' in the magic basis: 

cnot£ 2) = ^ 

\ — i 

(1 2) T (1 21 

and the eigenvalues and corresponding eigenvectors of (CNOT^^' CNOT M ' B ) are given by 



l>\ = S'2 : : L l«l) : — ! I , 1 02 ^ j 

(34) 

= M4 = -1, |a 3 ) = 75 I 1 1 |«4) = 
We now choose \v\) = [0,3), fa) = [0,4), fa) = fa) and fa) = fa), and calculate the matrix K2 from Eq. (fTUf : 

^ N ° T = -~! V n V n I <■>•-' 
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Next, we calculate the matrix A from Eqns. (fTTj) and ([T]). Choosing Ai = A2 = 1, A3 = A4 = i, we find the solution 
a = n/4, (3 = ^ — and hence 

/ 1 i 

A""-^ SllSl. (36) 



i 1 



Finally, we obtain iCi from Eq. (|l'2p: 



^ CNOT = o 1 _ • \ ; ; 1 ■ (37) 




We now decompose each of the matrices if^ 01 ", A CNOT and i^CNOT [Eq ns _ ([35]) - (|3T[) ] into a sequence of single-qubit 
rotations and (SWAP) Q -gates, using the procedure described in Section Til CI Starting with K§^ OT and comparing 
the right-hand sides of Eqns. ([35]) and (fT3|l for condition (|14a|l . we obtain: 

^CNOT = i? (l) ( | ) _ (38) 

Since yl CNOT — ^swBcii^ decompositions of these matrices are the same: 



a cnot = R (i) ^ ys^AP r£> (tt) VSWAP. (39) 
Finally, the decomposition of the matrix K™ OT £ SU(2) ® SU(2) is obtained using Eq. (HHJ), and we find: 

x CNOT = i? ( 1)( _| )i? ( 1)( | )il ,( 2 ) ( | ) , (40) 

The total decomposition of CNOT*- 1 ' 2 -' is then given by: 

CNOTg$ Ap) = ^CNOT^CNOT^CNOT 

= R^H~)R^Hl)R^H~)VSWAPR^\n)VMAPR^\l). (41) 

The number of operations in Eq. (|41[) cannot be further reduced by applying commutation relations or Euler's theorem. 
Including possible permutations of rotation angles, we obtain the general form of Eq. (|4"Tj^ : 

CNOTg$ Ap) « = i#)(a|^ (42) 

for a' , a") £ {(1, 2, —1, — a), (2, 1, a, 1)}, with a = ±1. Using the B-gate, the shortest-sequence implementation 
of the CNOT^'^-gate is given by: 

CNOT^ = i?«(a|)i?i 1 )(a'|)i?i 2 )(- ( xV|)S(a'|,a'|)i?W(-a|), (43) 

for a, a' = ±1. The sequences Eqns. (J3DJ), (EU, @2J) and |g3J) for implementing f7 st -* Bc11 and the CNOT-gate are 
building blocks for the implementation of quantum operations for three or more qubits^. 



C. The quantum Fourier transform 

The quantum Fourier transform (QFT) lies at the heart of Shor's factoring algorithm^! and is given by (for two 
qubits)2S: 

Viz 1 i I 
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So far, the QFT has been implemented in an NMR system^ 6 ., using ion qubits^ 7 . and using phononic qubitsS^, but not 
yet with solid-state qubits. J- 2 = CNOT^ 2,1 ) is a maximally entangling gate, but T itself is not*22. In order to find 
the shortest sequence of operations required to implement J 7 , we thus cannot use our optimization protocol for gates 
with maximally entangling capacity [as defined and described in Section [II Bj . but need to use the method for general 
gates U € SU(4) outlined in Section Hi Al To this end, we use the decomposition of T into the controlled phase gate 
CP plus two Hadamard gates proposed by Coppersmith et 



i 30. 



T = SWAP • H {2) ■ CP ■ H (1) 



with 



(45) 



= -4= 



CP 



H (i) 




(46a) 



(46b) 



(46c) 



The Hadamard gates ifW and can be decomposed into single-qubit rotations using Eq. (fTS"]) and we find 



ff« = R<i\-K)Rf{^) 



1,2. 



(47) 



We now proceed to calculate the decomposition K1AK2 of CP and to translate each of these matrices into single- 
qubit rotations and two-qubit interactions according to the protocol given in Sec. [XTJ Starting with CPmbi the 
representation of CP in the magic basis, we obtain from Eqns. (|46b|) , ([H]) and @: 



CPmb — 7j 



1+i 1 + i 

2 

2 

1 + i 



The eigenvalues and eigenvectors of CP^ ib CPmb are given by 



Mi = M2 



fJ-3 = M4 







f 









(!) 


. I fl 2) = 









Vi 


(V 




1 




j l«4> = | 










\0, 




V0 



(48) 



By choosing \v% 



for i = 1, 



. 4 and using Eq. (fT0"|) , we obtain the matrix K§ p : 
/I 

K? p = - 




(49) 
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where we have used Eq. (fT5|) to obtain the decomposition of K^ p ■ Assigning Ai = A2 = \fi and A3 = A4 = 1, we 
obtain from Eq. (fTTjl the solution a = ir/8 and /? = 7 = 0, and thus from Eq. ((TJ: 



(50) 



cos(|) isin(f) 

,CP _ I cos(f) isin(f) 

isin(f) cos(|) 

isin(g) cos(f) 

= i?i 1) (7r)(SWAP)3Ri 1 )( 7 r)(SWAP)3. (51) 



The decomposition (|5ip of A CP is obtained by calculating the concurrence (see Eq. Q61[l) of the column vectors of 
A CP , which yields C = 4cos 2 (7r/8) sin 2 (7r/8) = 1/2. Following the method outlined in Sec. Ill Cl we then calculate the 
corresponding time t* such that (SWAP) Q (* ' yields the same value 1/2 of the concurrence (amount of entanglement) 
when applied to the standard basis vectors. We find that two applications of the (SWAP) ^-interaction, each followed 
by a rotation of the first qubit over it around the x-axis (which interchanges the first and the third row, and the 
second and the fourth row of the matrix to which it is applied) yield a decomposition of the matrix (|50[) . 
Having found the matrix A CP , the matrix K^ p and its decomposition can be calculated from Eqns. (|12p and (| 15[) . 
respectively, 



Ki P 



/111 1 \ 

1 ■ 1 - 1 ■ 1 ■ 

\ i —i—i i j 

= Rf\^)RfX-l)R { ?{-\)R£\-l) 

= R^i~)R^(~)R^(~)R^( 3 iy (52) 

The decomposition of the quantum Fourier transform T is then given by: 

T = SWAP ^) ( -J)i ? W(_|)4i)(^) ( SWAP)ii?W( 7 r)(SWAP)ii?( 1 )(-|)i?( 2 )(-|^ 

= RW { ^)R(i)(-l)R(2) (^) ( SWAP)fi?( 1 )( 7 r)(SWAP)ii?W(~|)i?( 2 )(-|)i?( 2 )(|). (53) 

Using commutation relations to minimize the number of rotations in Eq. (|53p , we obtain the decomposition of T that 
contains the smallest number of (SWAP) Q - and single-qubit gates: 

^(swap)= = Rf{~\ )i?f(J) (SWAP)f R^(n) (SWAP)i R^{^)R^{^), (54) 



and analogously we find for the optimal decomposition of T using the B-gate^L: 

^b = R^H~)b(j,^)R^K^)R^H^)B(~,^)R^H^)R^(^)R^(^)- (55) 



■^(SWAP)° and Tb also give alternative implementations of the CNOT^ 2 ' 1 ' gate, which consist of applying the sequences 
([54")) or ([55]) twice. These, however, are longer than the sequences obtained by optimizing CNOT directly [Eqns. (|4"2"]) 
and (|43p ]. since they do not implement the transformation from standard to maximally entangled basis with the least 
amount of two-qubit operations. 



IV. APPLICATION TO ELECTRON SPIN QUBITS AND SUPERCONDUCTING FLUX QUBITS 

For solid-state qubits, implementation of two-qubit quantum gates in an efficient way is important because of 
their short coherence times. In this section we discuss the implementation of the above results for electron spin 
qubits and superconducting flux qubits and estimate relevant time scales. An electron spin qubit consists of a single 
electron confined in a quantum dot, a small island in a semiconductor structure that can be filled with electrons in 
a controlled way22. The qubit is encoded in the spin degree of freedom. Single qubit rotations can be generated 
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by applying an electron spin resonance (ESR) pulse to the electron which is described by the evolution operator 
Un{t) = exp[— (i/H) J^Tifi(T)dT], with Ttn(t) — —^fvyB(t) ■ a. Rotations of qubit 1 and 2 are then represented as: 

i#>(/9) = e -* ifMt -* 91 and R.W(0) = e -*</w®** (56) 

with p{t) = —7 J B(t)(1t, 7 the gyromagnetic ratio, B(t) the magnetic ESR field (applied in the direction perpen- 
dicular to the Zeeman-splitting field), n = (sin 8 cos (j>, sin sin 0, cos ff) a unit vector on the Bloch sphere [0G[0, 7r), 
</>€[0, 27r)], a = (cr x , a y , a z ) and / the identity matrix. Two-qubit interactions arise from the Heisenberg exchange 
interaction, so that the (SWAP) a gate [Eq. (fT7|) ] is the natural two-qubit gate for electron spin qubits. Both single- 
spin rotations and v SWAP-operations have recently been demonstrated for spin qubit a 33 ' 34 , and typical times are 
W ~ 100 ns (for magnetic fields of 1 mT) and t^j^ ~ 180 ps. 

Flux qubits^ consist of a superconducting loop interrupted by three or four Josephson junctions. The qubit 
basis states consist of the direction of the current that is circulating around the loop and single qubit rotations are 
generated by applying resonant microwave radiation. These rotations are described by the same expression (|56p as for 

spin qubits with {3 — > jr L Qj(r) cos 4>j(T)dT for a x (rotations around the x-axis) and (3 — > — ^ J fij(r) sin (j)j(r)dT 
for cjy (rotations around the y-axis), j = 1,2. Here f2j(r) and 4>j{r) denote the amplitude and phase of the applied 
microwave signal. Measured Rabi oscillations arc of the order of 1-10 ns 3e . Various proposals to achieve a tunable 
coupling mechanism for flux qubits have recently been put forwar d 19 ! 37 ' 38 and/or realized 39 . In each of these the 
proposed two-qubit gates are equivalent to the B-gate^. In Refji 9 - the creation of controllable coupling between two 
detuned flux qubits via the quantum inductance of a third flux qubit is suggested, as described by the Hamiltonian fll8p . 
Predicted B-gate operation times for this system range from ~ 10 ns 1 ^ down to ~ 2 ns^ 7 -. 

In practice, manipulating qubits introduces decoherence. The single-qubit decoherence time T2 has not been 
measured yet for electron spins, but recent experiments show that single-spin Rabi oscillations remain visible for up 
to 1 /is (where each oscillation takes ~ 100 ns)^ 3 - and ensemble decoherence times T 2 * > 1/is 34 , so that T2 is expected 
to be > l^s. For flux qubits, measurements of T2 range from 15 ns 3 ^ to a few /j,s— , while a Rabi oscillation requires 
~ 5 ns. T2 is thus sufficiently long to observe at least a few Rabi oscillations. It is unknown, however, to what degree 
decoherence will affect a quantum gate operation which consists of several single qubit rotations on different qubits 
as well as two-qubit interactions: minimizing the total number of operations used 4 ^ is then likely to be an essential 
factor for achieving high gate fidelities. 

V. CONCLUSION 

We have presented a systematic protocol for calculating the optimal (defined as consisting of the smallest number 
of single- and two-qubit operations) implementation of a desired two-qubit gate U S SU(4) in terms of a given 
tunable two-qubit interaction Ui n t(t) and single-qubit rotations. We have illustrated the decomposition method by 
calculating the shortest sequence of operations required to implement the transformation from the standard basis to 
the Bell basis, the CNOT gate, and the quantum Fourier transform, using either Heisenberg exchange interaction 
or quantum inductive coupling. The general method presented here is a useful tool to find the smallest number of 
operations that are needed to realize a desired two-qubit gate for any type of qubit, using the single- and tunable 
two-qubit operations that are naturally available in the qubit system. 
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VII. APPENDIX 



In this appendix we calculate the concurrence of the column vectors of the matrix 

t A(e)^ i )( M )i4 2 )(^)4 i )(e)4 2 )(o, 



(57) 
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which is used in the proof of Proposition 1 in Sec. Ill Bl By writing out Eq. (|57f explicitly, we obtain 



T 




I B X C X + B 2 C 2 - iB* 2 C 3 + iB\ C* 4 
a(B 2 d + BxC 2 + ^B^C , 3 - iB* 2 C A ) 
a{-iB*d + iB\C 2 + B^ + B 2 C 4 ) 
V iBld - iB*C 2 + B 2 C 3 + B X C 4 



(58) 



with 



B 2 
Cx 
C 2 
C 3 

c 4 



[I v 

cos — cos — 
2 2 

. P v 
sin — cos — 
2 2 

cos % COS ~ 
2 2 

i . c 

cos- sin - 

. t c 

sin -cos - 

. e . c 

Sin -Sin - 
■ 7r ^ 



/i i/ 
■ « sin — sin — 
2 2 
/i . v 
i cos — sin — 
2 2 



exp(- 



(59) 



The concurrence C of each of the column vectors \wi) of T is found by calculating the square roots Ai > A2 > A3 > A4 
of the eigenvalues of the matrix 



with pi = \wi)(wi\, i = 1, 



M = Pl (<Jy (g) <Ty)p* {<Jy <g) Oy ) , 

, , 4, and substituting 

C = max{Ai — \ 2 — A3 — A4, 0}. 



(60) 



(61) 



The conditions on e, fx, v, £ and £ in Eq. (I14[) are then obtained by evaluating C = 1. The resulting expressions are 
in general lengthy and therefore not given here. As an example, consider the special case £=C=0. In this case, T 
[Eq. (pTT))] reduces to 



/ 



T, 



?=C=o 



Bi 
aB 2 



-iaB 2 %aB\ 



B 2 -iB* iB\ 
aB\ iaB\ —ictB^ 
aB\ olB 2 



(62) 



V %B\ 



-iB% B 2 



We calculate the concurrence of the last column vector \1ju4) of Eq. (f62|) (the same result is obtained for the other 
column vectors). Let 



p = |w;4) (w 4 1 



Then M [Eq. (JBD])] is given by 



-uB x B* 2 
—iaB\B 2 
-%B\ 



-a*B\B 2 ia*B\B* 2 i{B\f 



iu*B x B 2 



a*BxB* 



-iaB\B* 2 
a.B\B 2 



(63) 



\B Y \ 2 P -B{B 2 Q iB'{B 2 Q i{B{) 2 P 

1 I ~B X B* 2 Q \B 2 \ 2 R -i(B*) 2 R —iB*B%Q 

2 I -iBxB 2 Q iB 2 R \B 2 \ 2 R B\B 2 Q 
~iB 2 P iB x B 2 Q BxB*Q \Bx\ 2 P 



(64) 
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with 

P = \B X \ 2 + (a*f\B 2 \ 2 

Q = alBil 2 + a*\B 2 \ 2 

R = a 2 |Bi| 2 + \B 2 \ 2 . 

The eigenvalues of Eq. (f64|) are the solutions of the equation 

A 3 (A-2(|S 1 | 2 P+|S 2 | 2 i?)) =0 
^A = or A = iBil 4 + 2cos(e7r)| J Bi| 2 | J B 2 | 2 + \B 2 \ 4 > (65) 

so that 



C = v /|S 1 | 4 + 2cos(e^)|i? 1 |2|B 2 |2 + |i ?2 |4. ( 66 ) 

Finally, 

C=1 {^0 and (|Bi| a = 0,|fl a | 2 = l or l^ 2 = 1, \B 2 \ 2 = 0). (67) 

The latter solution corresponds to fx, v S {0, 7T, 27r, 37r, . . .}. Eq. (|6"T)) corresponds to the solution given in Eq. (fT4")) . 
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